Energy transduction in periodically driven non-Hermitian systems 
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We show a new mechanism to extract energy from non-equilibrium fluctuations typical of peri- 
odically driven non-Hermitian systems. The transduction of energy between the driving force and 
the system is revealed by an anomalous behavior of the susceptibility, leading to a diminution of 
the dissipated power and consequently to an improvement of the transport properties. The general 
framework is illustrated by the analysis of some relevant cases. 
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In the last years a growing interest in mechanisms 
for energy transduction by rectification of unbiased ther- 
mal fluctuations has arose, partly motivated by problems 
from cell biology. Several phenomenological models, cat- 
egorized as thermal ratchets or Brownian motors, have 
been proposed. These engines operate at molecular level 
[H-H , although their potential implementation in a larger 
scale would be of evident interest. In this context, a num- 
ber of methods for particle separation have been recently 
proposed based on several variants of the ratchet concept 
|3|,|6|. Motivated by this interest, we investigate how to 
take advantage of non-equilibrium fluctuations to opti- 
mize the energy consumption. We propose a new mecha- 
nism for energy transduction in a class of non-equilibrium 
systems when they are acted upon by a weak periodic 
force. The coupling between the external driving and 
the out-of-equilibrium fluctuations leads to a minimiza- 
tion of the dissipated power. The minimum occurs when 
the frequency of the external force matches a characteris- 
tic frequency of the system, thus manifesting a resonant 
behavior. 

Let us consider the class of differential equations 



fit*(2,t) = (£o + A(i)£i)*(f,t), 



(1) 



describing the dynamics of a probability density or of a 
hydrodynamic field ^(x, t), where x represents a coordi- 
nate. The dynamics is governed by the non-Hermitian 
operator Cq and is influenced by the action of a peri- 
odic force which introduces the perturbation X(t)C\, with 
\{t) = X e tuJt . 

We will analyze the response of the system to the ex- 
ternal perturbation by using the standard procedure of 
the linear response theory. Accordingly Eq. (nl) can be 
solved formally yielding 

$>(x,t) = e (t - to)c °y Q (x) 

+ f dTA(T)e (t ~ T)£o (£ie (T - to)£o * (x)) (2) 
Jt„ 

= tf(°>(z,t)+Atf (£,*), 

where ^q{x) is the initial condition, which corresponds 
to the stationary state of Eq. (pi) in the absence of the 
external force [0J|]. Expansion of the term Cie tCoi &o(x) 



in series of the eigenfunctions of the operator Cq, 4>n with 
eigenvalue a n ; n = 0,1, ... 

oo 

C ie tc ^ (x) = ]>>„</>„(£) + «(£)}, (3) 



where c n are the corresponding coefficients in this expan- 
sion, then leads to 



Ax{t) = / xA^(x,t)dx = / drx(t - r)A(r), (4) 

which defines the susceptibility x(£). 

We will assume the existence of a dominant time scale 
governing the relaxation process, which corresponds to 
the n = 1 mode in the expansion (0) . Since the remaining 
modes decay faster we can truncate the series retaining 
only the first term. Thus, considering only contributions 
of the first mode, 



X(t) = Ae ait +A*e<\ 
with A defined as 

A = c\ I x<pi(x)dx, 



(5) 



(6) 



the explicit expression of the susceptibility as a function 
of the frequency follows from the Fourier transform of 
Eq. (pi) 



xM = 



A 



1 



,4* 



1 



h(3- i(a + 1) h /3 - i{a - 1) 



(7) 



where a\ = R\ + Hi, j3 = Ri/Ii and a = lu/Ii, with R\ 
and I± being the real and imaginary parts of ai, respec- 
tively 

In Fig. 1, we have plotted the modulus of the sus- 
ceptibility as a function of a for different values of [3. 
During the relaxation process of non-equilibrium fluctu- 
ations, the susceptibility undergoes a resonant behavior 
when the frequency of the force matches the imaginary 
part of the first eigenvalue of the non-perturbed opera- 
tor Cq. This behavior reveals the resonant coupling be- 
tween the periodic force and the non-equilibrium source, 
responsible for the non-Hermitian nature of Co ■ 



The appearance of this resonance leads to a diminution 
of the dissipation in the relaxation process of the fluctu- 
ations of x. To illustrate our assertion, let us suppose 
that Eq. ([!]) represents the non-Hermitian Fokker-Planck 
equation 



P(w) 



2tt 



2-k/u 



dtp,, 



(15) 



i) 



For the particular case of the Brownian particle, Pf 
reads 



d t ^(x, t) = -Va • (v{x)^ - DV S V - b\(t)VV s U(x)), (8) Pf = p (t) 

where v(x) is a non-potential drift, b a mobility, D = 
KsTb the corresponding diffusion coefficient and U(x) 
the potential related to the external force. In this case 
the operators Cq and L\ act on the field as 



d(x) 
~dT 



J2WRe{F n (t)-J*}. (16) 



n^O 



£o* = -V a • (v(x)V - DV^), 
£i* = Vs ■ (bW s U{x)). 



In Fig. 2a we have plotted the corresponding P(w). To 
this end we have assumed that the motion of the particle 
takes place in two dimensions, with periodic boundary 
conditions, thus, the set of eigenfunctions 4>t{x) are the 



(9) 



Fourier modes with eigenvalues as = —Dk 2 



' ■ k. We 



Amongthe physical realizations of the model described 
by Eq. (0) we can quote the case of a Brownian parti- 
cle advected by a constant drift v acted upon by a force 
F(x*, t), or a field-responsive particle in a vortex flow un- 
der the influence of an oscillating magnetic field fipj . 

We are interested in analyzing the energy dissipated 
by the system in the dynamic process governed by Eq. 
(0). The dissipated power is 



dxJ ■ Vj/j, 



(10) 



have considered a force of the form F(x,t) = A(£)(l + 
cos(fc • x))v, where v is the unit vector pointing along 
the direction of the drift. The figure shows that P(ui) 
achieves its minimum value at the resonant frequency. 
The negative character of this quantity indicates that 
the system is acting as a generator. 

In Fig. 2b we have represented that quantity for the 
more complex case of the mesoscopic dynamics of a field- 
responsive Brownian particle in a vortex flow with vortic- 
ity wo, under an oscillating magnetic field, H(t), P,fl0|. 
In this case, the power supplied by the external force is 



where J = -£>Vj* - b\(t)^y S U [x) is the diffusion cur- 
rent and \i = KBT\n^(x, t) + \(t)U(x) the correspond- 
ing chemical potential. It splits up into two contribu- 
tions, one coming from the underlying diffusion process 
and another accounting for the power supplied by the ex- 
ternal force. Since we are interested in the effects of the 
external field, we focus on the externally supplied power, 

P F = -X(t) / dxJ ■ V S U = / dxJ ■ F(x,t). (11) 

By expanding F(x, t) in series of the eigenfunctions of Cq 
F(x,t)=F Q (t) + J2(Fn(t)(f>n(S)+F:(t)cf>* n m, (12) 

and substituting this equation into Eq. (O) we achieve 

'd(x) 



Pf = H{t) 



fd(M) 



dt 



- uj x (M) 



(17) 



P F =F (t) 



dt 



(v(x))\+Y,WRe{F n (t)-j:}, (13) 



where M is the magnetization. This figure shows that in 
this case the transduction of energy occurs in two differ- 
ent regimes. In the low frequency regime Pq, the power 
dissipated in the Debye relaxation [O] , achieves negative 
values while P c , the energy dissipated due to the coupling 
between the drift and the external force, takes negative 
values for high frequencies. In both cases the averaged 
dissipated power P(uj) exhibits a minimum value at the 
resonant frequency, showing the resonant character of the 
mechanism for energy transduction. 

The analysis of the energy dissipation allows us to 
study the transport properties. The presence of an ex- 
ternal driving, forces the system to move with a mean 
velocity, v m , different from the drift v(x), thus leading to 

the appearance of a drag force, Fd = — %■ v m [ fL3] , with k 
a friction tensor accounting for dissipation in the system. 
The resulting dissipated power is 



where fl is the volume of the system. To obtain Eq. <Rc 
we have used the result 



J — Vm. ' K ' Vm,- 



(18) 



dx J 



d(x) 
dt 



(vm- 



(14) 



In the case of the Brownian particle advected by a con- 
stant drift, the friction tensor is a scalar. Its expression 
follows form Eq. (§X 



This expression follows from the Fokker-Planck Eq. (g) 
through the definition of J. 

The quantity of interest is the dissipated power aver- 
aged over the period of the external force, 



1 



k b 



-P, 



(19) 



whose behavior is essentially shown in Fig. 2a. 



For the field-responsive particle the dissipated power 
given through Eq. (|17|) consists of two independent 
contributions corresponding to longitudinal (Pq) and 
transversal (P c ) effects, with respect to the direction of 
the magnetic field. Consequently, associated to the last 
one, which corresponds to the viscous dissipation occur- 
ring when the magnetic field acts on the fluid, we can 
define a friction coefficient as 



i 



kf 



T\2 



(vl) 



uj a • «M> x H). 



(20) 



where v^ = cuga is the transversal component of v m , with 
a the radius of the particle and loq the modulus of the 
vorticity. 

It is interesting to analyze the behavior of this quan- 
tity in terms of the parameter /3, which is given in this 
case by D r /coo, with D r , the rotational diffusion coeffi- 
cient. As can be seen in Fig. 3, as (3 grows the friction 
coefficient kf becomes positive in the entire frequency 
range. In fact the frequency for which up = goes to in- 
finity in the limit of a Hermitian dynamics. Fig. 1 shows 
something analogous, that is, the resonance disappears 
when [3 grows. Thus, we conclude that the phenomenon 
described is genuine of non-Hcrmitian systems. 

The anomalous behavior exhibited by the friction co- 
efficient is a direct consequence of the form of the suscep- 
tibility. In equilibrium, the fluctuation-dissipation theo- 
rem, implying that uj\mxx(u) > 0, manifests that during 
the relaxation of the fluctuations around an equilibrium 
state the system always dissipates energy. Nonetheless, 
in the non-equilibrium case shown in Fig. 4 the imag- 
inary part of the susceptibility achieves negative values 
for positive frequencies thus violating the aforementioned 
inequality. In this figure, obtained for the particular case 
of the Brownian particle in a constant drift, the analyti- 
cal results are compared with numerical results from the 
corresponding Langevin equation by means of a second 
order Runge-Kutta method [ff4j . For the field-responsive 
particle, we obtain the same behavior, which essentially 
corresponds to Po, plotted in Fig. 2b. This fact indi- 
cates that the system is generating energy instead of dis- 
sipating power, which manifests at a macroscopic level 
through the diminution of the friction coefficient. An 
anomalous behavior of the response was first discussed 
in the context of current generation by noise induced 
symmetry-breaking in coupled Brownian motors pjfl . In 
addition, several similar phenomena have been reported 
in |lqjl7|l . These papers were focused on specific models, 
raising the question under which circumstances this be- 
havior arises. In this letter, we have found the conditions 
for these phenomena to occur in a quite general class of 
systems. 

Notice that our results differ from the ones obtained 
when LRT is applied to non-equilibrium Hermitian sys- 
tems. For example, in Ref. M fluctuation-dissipation- 
type relationships are derived, by assuming that the re- 
laxation occurs as in the case of an equilibrium state. In 



the present context, this assumption does not hold, since 
we are dealing with non-Hcrmitian systems whose eigen- 
value spectrum is complex. Consequently, perturbations 
do not relax exponentially as they do in equilibrium. 

The theoretical framework discussed in this paper can 
be applied to a number of problems formulated in terms 
of non-Hermitian dynamics. Among them, one could 
mention the transport of classical particles advected by 
a quenched p8| velocity field. This process, governed by 
a Fokker-Planck equation with random drift, models the 
diffusion in porous media. It has been recently shown pq] 
that, when the velocity field displays correlations in both 
longitudinal and transversal directions, the eigenvalues 
occupy a finite area in the complex plane. Consequently, 
this system evolves according to a non-Hermitian dy- 
namics. If the particles can respond to an external field 
the same phenomenology described in this paper holds, 
i.e., as a consequence of the diminution of the dissipa- 
tion transport through the porous medium becomes en- 
hanced. The phenomenon we study may also arise in 
population biology problems, in particular in the gener- 
alization of the Malthus-Verhulst growth model proposed 
by Nelson and Schnerb p9| . The linearization of this 
model around its steady state yields a non-Hermitian evo- 
lution equation. When a periodic driving is introduced 
as a time dependence of the resources of the medium, 
the minimum achieved by the dissipated energy is now 
related to a resonant optimization of these resources. 

Additionally, our approach might unify the explana- 
tion of other resonant transport phenomena previously 
reported. The Senftleben-Beenakker effect observed in 
gases of polyatomic molecules shares the phenomenology 
inherent to our model. This effect occurs when the gas 
is under the action of both a constant magnetic field and 
an oscillating field parallel to the first one pCfl . Larmor 
precession causes the non-Hermiticity of the Boltzmann 
equation describing the dynamics of this system. The 
resonant frequency is related to Larmor's frequency. Un- 
der these conditions, the viscosity of the gas manifests a 
non-monotonous behavior as a function of the frequency 
similar to the one depicted in Fig. 3. Another exam- 
ple exhibiting analogous characteristics is the negative 
viscosity effect observed in field-responsive fluids |)|,|l()f , 
under a non-potential flow and submitted to an AC field. 
This effect consists of a diminution of the viscosity due 
to the presence of the periodic field. The field-responsive 
phase acts as a transmitter of energy between the exter- 
nal force and the system, leading to an improvement of 
the transport. 

In summary, we have proposed a mechanism for en- 
ergy transduction in non-equilibrium systems, based on 
the possibility of extracting energy from the relaxation 
process of out-of-equilibrium fluctuations. Under the ac- 
tion of an oscillating force, systems which evolve accord- 
ing to non-Hermitian dynamics act as transducers. Con- 
sequently, the energy dissipated in the system dimin- 
ishes achieving its minimum value when the frequency 
of the external driving matches the resonant frequency. 



This diminution of the amount of energy dissipated has 
a strong influence on the macroscopic properties of the 
system, leading to an enhancement of the transport or 
more generally to an optimization of the consumption 
of energy. Due to the intrinsic non-equilibrium nature 
of the fluctuations, energy transduction does not require 
any further ingredient, as occurs in ratchet-like engines, 
in which the presence of a parity-symmetry-breaking po- 
tential is an unavoidable condition for transferring en- 
ergy. 
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FIG. 1. Non-dimensional modulus of the susceptibility as 
a function of the parameter a. Continuous line corresponds 
to j3 — 0.1, dotted line to /3 — 0.5 and dashed line to (3 — 1. 
The resonance fades away practically for f3 « 10. 
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FIG. 2. 2a.-Non-dimensional dissipated power for a parti- 
cle under a constant drift as a function of the non-dimensional 
parameter a.2b.-Same for a field-responsive fluid. The dotted 
line corresponds to the total dissipation -P(w). 
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FIG. 3. Non-dimensional friction coefficient as a function 
of the parameter a for different values of (3. Continuous line 
corresponds to j3 = 0.1, dotted line to /3 — 0.5 and dashed 
line to P — 1. 




FIG. 4. Non-dimensional imaginary part of the susceptibil- 
ity for a Brownian particle advected by a constant drift as a 
function of the parameter a. Solid line corresponds to the an- 
alytical computation whereas dots have been obtained from 
numerical simulations. 



